Ultracold heteronuclear molecules and ferroelectric superfluids 
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We analyze the possibility of a ferroelectric transition in heteronuclear molecules consisting of 
Bose-Bose, Bose-Fermi or Fermi-Fermi atom pairs. This transition is characterized by the appear- 
ance of a spontaneous electric polarization below a critical temperature. We discuss the existence of 
a ferroelectric Fermi liquid phase for Fermi molecules and the existence of a ferroelectric superfluid 
phase for Bose molecules characterized by the coexistence of ferroelectric and superfluid orders. 
Lastly, we propose an experiment to detect ferroelectric correlations through the observation of 
coherent dipole radiation pulses during time of flight. 

PACS numbers: 03.75.Kk, 03.75.Ss, 32.10.Dk, 33.15.-e 



Arguably one of the next frontiers in ultracold atomic 
and molecular plwsics is the study of two-species atomic 
mixtures P, S S IE 1 ^'^^ ultracold heteronuclear 
molecules composed of two-species alkali atoms such as 
KRb 0,1113, RbCs [13 and NaCs [ij. This frontier 
may be advanced through the use of Feshbach resonances 
which have already been observed in mixtures of two- 
species alkali atoms 0, 0, Q, and may also become 
a crucial tool for tuning physical properties of heteronu- 
clear systems. 

Ultracold heteronuclear molecules made of Bose-Bose, 
Bose-Fermi or Fermi-Fermi atom pairs offer many new 
opportunities compared to standard (Bose or Fermi) 
atomic systems because of their additional degrees of 
freedom. For instance, when heteronuclear diatomic 
molecules are formed from neutral atoms, electric charge 
is transferred from one atom to the other leading to an 
electric dipole moment |p| = Q{d)d, where d is the sepa- 
ration and Q{d) is the effective charge transfer between 
constituent atoms 15]. These electric dipoles have equal 
magnitudes but random orientations at high tempera- 
tures leading to a vanishing average electric polarization. 
However, at low temperatures, the dipoles may all point 
to a particular direction producing a finite average elec- 
tric polarization density (P), characteristic of a ferro- 
electric state. In addtion, when ultracold heteronuclear 
molecules form a Bose-Einstein condensate (BEC) a fer- 
roelectric superfluid state proposed in this manuscript 
may be accessible experimentally. 

Our main results are as follows. When heteronuclear 
molecules are composite fermions (Bose-Fermi pairs), a 
phase transition occurs separating a paraelectric Fermi 
liquid (PEL) from a ferroelectric Fermi liquid (EEL), as 
shown in Fig. ^a). These phases do not exist in stan- 
dard condensed matter systems, since all experimentally 
known ferroelectrics are very good insulators |l6l |. and 
are not Fermi liquids. Furthermore, when heteronuclear 
molecules are composite bosons (Bose-Bose or Fermi- 
Fermi pairs), a ferroelectric transition with critical tem- 
perature TpE may occur either above or below the BEC 
temperature Tbec- When the molecular dipole moment 
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FIG. 1: Schematic phase diagram for polar a) Fermi, and b), 
c) Bose molecules, where (P) and {(j)) are the corresponding 
order parameters for ferroelectric order at Tfe, and for BEC 
at Tbec, respectively. 



and/or density are sufficiently large then TpE > Tbec 
leading to a paraelectric Bose liquid (PBL) for T > Tps, 
a ferroelectric Bose hquid (FBL) for TpE > T > Tbec, 
and to a ferroelectric BEC (FBEC) for T < Tbec, a-s 
shown in Fig. ^b) . When the molecular dipole moment 
and/or density are sufficiently low then Tps < Tbec 
leading to a paraelectric Bose liquid (PBL) for T > 
Tbec, a paraelectric BEC for Tbec > T > TpE, and 
to a ferroelectric BEC (FBEC) for T <Tfe, as shown in 
Fig.^c). The FBEC phase corresponds to a ferroelectric 
superfluid. 

Hamiltonian: To describe all these phases of ultracold 
heteronuclear molecules, we use the Hamiltonian density 



H{r) = 



2M 



^i + V(r) 



+ ^E/ dr'n^,{r)Ud{r,r'M) 

-\T.I rfr'P„(r')Jo,/3(r,r')P/3(r), (1) 
describing a weakly interacting gas of dilute po- 
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lar molecules, where creates a Bose (Fermi) 

molecule at position r with pseudo-spin-cr, and na-{r) = 
tjj*{r)tpa{i^) is the density operator. The first term in 
Eq. ^ is the kinetic energy, where ^ is the chemical 
potential, and ^(r) is the trapping potential. The sec- 
ond term is the density-density interaction between Bose 
(Fermi) molecules, and for a contact interaction is given 
by [/c;(r,r') = Ud{Ar), where U — Ana/M is the ampli- 
tude of the interaction and (5(Ar) is the delta function 
with Ar = r — r'. Here, a is the scattering length of 
the corresponding Bose (Fermi) molecules, and M is the 
molecular mass. The third term is the electric dipole- 
dipole interaction between molecular dipoles at positions 
r and r', where P(r) — p(r) '^^(r) is the polarization 
operator, and p(r) is the molecular dipole moment at 
position r. The electric dipole-dipole interaction is given 
by J(r,r') = X;a^/3 P«(r) Ja,,3(r, r')p/3(r'), where 



>^a,/3(r,r') 



1 3[Ar„Ar;3]/|Ar|^ 



0a,l3 



47reo 



lArP 
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Here, Sa^p is the Kronecker-delta and {q!,/3} ~ {1,2,3} 
label the vector components. 

Electric polarization: Within the Hartree-Fock approx- 
imation, the Hamiltonian density reduces to 



M^) - ~—~ll + V{v) + Un,{v) 

- J2[ *'p,(r)J„,^(r,r')Po,/3(r'), (3) 

where the electric polarization density Po(r) = 
p(r)no(r), with no(r) = X)<t "'0,o-(i") being the local den- 
sity of Bose (Fermi) molecules at r. The spatially aver- 
aged polarization (P) = J dr'Po{r) /Vc can be rewritten 
as 



(P> = -^/rfrp(r)^|0,(r)|V.(eO, 



(4) 



where Vc is the sample volume, 0i(r) and are eigen- 
functions and eigenvalues of the Hamiltonian Ho{r), and 
fvi^i) — ^/i^^'^' — il] is the Bose (Fermi) distribution for 
Bose (Fermi) molecules when ry = 1 (—1). The solution 
of Eq. I0J is non-trivial, however analytical insight can 
be gained for homogenous systems where T^(r) =0 and 
Po(r) = Pq is independent of r. We discuss the homoge- 
neous case first and then analyze the case of a harmonic 
trap. 

Ferroelectric critical temperature: In the ferroelectric 
state, all molecular dipoles are pointing along the same 
direction m such that p(r) = p = |p|m and Pq = |Po|ni- 
The critical temperature TpE for the ferroelectric tran- 
sition is found from the slope of Eq. Q with respect to 
Pq evaluated at Pq = 0, leading to 



1 + J(q 



dN 



0, 



(5) 



from which TpE can be calculated. Here, J(q) with q = 
k — k' is the Fourier transform of J(Ar), and N is the 
number of Bose or Fermi molecules. Introducing the Kac 
parameter 7 as a cutoff for short distances (|Ar| < I/7) 
leads to J(q) = X^a,/^ Pa'^a,/3(k, k')p/3, where 



J„^;3(k,k') = 



S{q) f 3qQq/3 



£0 



(6) 



Here, S{q) = sinq/q^ — cosq/q^ where q = |q|/7. Thus, 
we obtain an implicit relation for the TpE 



M 

3eo 



(7) 



in terms of the molecular density n — N/Vc and the 
isothermal compressibility k{T) — {l/n'^){dn/d^)T- 

The ferroelectric instability is accompanied by a diver- 
gence of the dielectric function in the long-wavelength 
and low- frequency limit. Using linear response the- 
ory, the dielectric function e(q, zw„) can be related 
to the density-density correlation function C(q, t) = 
{Trn{q, r)n(— q, 0)) and to J(q) via 



e(q, iwn ) 



Vc 



I ' dre™"''C(q,r), 



(8) 



where n(q, r) is the density operator. In the long- 
wavelength and low-frequency limit, C(q, t) is directly 
related to k{T) via the compressibility sum rule |l7| . 
Therefore, a divergent dielectric function occurs when 
e(q — ^ 0,iw„ 0) = leading to 1 + J(q 
0){dN/dn)T = 0, which is identical to Eq. ©. This 
relation can be applied to both Bose and Fermi systems. 
Next, we discuss TpE for a weakly interacting (a <g; A) 
and dilute {na^ ^ 1) gas of Bose and Fermi molecules, 
where a is the scattering and A is the thermal length. 

Fermi molecules: As a first application of Eq. j?)), 
we analyze Tfe for a weakly interacting gas of Fermi 
molecules at any T. In this case, the molecular density 
is given by n = 2F3/2{z)[l ~ 2Fi/2{z)a/ Xp]/ X% where 
\f = [27r/(Mr)]i/2 is the thermal length, and lead- 
ing to k{T) « [Fi/2(z)/F3/2(^) - 2F_,/2iz)a/XF]/{nT) 
where < z = exp(/3/i) < 00 is the fugacity, and 
F^(z) = [l/r(i/)] x^-^dx/lz-^e^ + 1] is the Fermi in- 
tegral. Here, r(i/) is the Gamma function. Thus, we 
obtain 



Tfe 
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Sen 
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^3/2 (^c) 



Ar 



n, 



(9) 



where Zc = exp{ii/TFE)- Notice that, in the classical 
{zc < 1) limit, Eq. © reduces to Tfe = 2|pp(l - 
2zca/ XF)n/ {3eo), which shows that Tfe oc n for a 
gas of classical electric dipoles. The Tfe for an ideal 
(non-interacting) gas of Fermi molecules can be ob- 
tained by setting a — 0. For an ideal gas, when Tfe 
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is much smaller than the Fermi energy ep, we obtain 
Tfe ~ (2V^/Tr)eF[\p\^n/{soeF) - which is valid 

for |p|2n/(eoef ) > 1 and [\p\hi/ (eoep) ~ 1]^/^ < 1. For 
T > Tfe a PFL phase exists and tor T < Tfe a FFL 
phase appears as shown in Fig.QJa). 

Bose molecules for T > Tbec- As a second ap- 
plication of Eq. Q, we analyze Tfe for a weakly in- 
teracting gas of Bose molecules when T > Tbec- 
In this case, the molecular density is given by 
= B3/2{z)[l - 4Bi/2(z)a/AB]/A| where As = 
[l/(27rMT)]^/^ is the thermal length, and leading to 
k{T) « [Si/2(z)/B3/2(z) -4B_i/2(z)a/AB]/(nT) where 
< z = exp(/3/i) < 1 is the fugacity, and By{z) — 
\l/T{v)]^^ x^^^dx/lz-^e^ - 1] is the Bose integral. 
Thus, we obtain 



T 



FE 



2|PP 

3£o 



Bi/2{zc) _ 4a 

B3/2{Zc) As 



S-i/2(^;c) 



n. (10) 



Notice that, in the classical {zc <C 1) limit, Eq. H1U|) 
reduces to Tfe = 2|pp(l — Azco/ \B)n/ [iso), which 
again leads to Tfe oc n for a gas of classical electric 
dipoles. The Tfe for an ideal (non-interacting) gas of 
Bose molecules can be obtained by setting a = 0. Notice 
that, the second terms in Eqs. © and H1U|) are different 
by a factor of 2 due to the degeneracy of pseudo-spin-1/2 
fermions in contrast to pseudo-spin-0 bosons. 

Bose molecules for T < Tbec- As a third application 
of Eq. |(7J), we analyze Tfe for an ideal (non-interacting) 
gas of Bose molecules when T < Tbec- In this case, 
the molecular density is given by n = _B3/2(^)/A'g -t- ris 
where Us = z/[Vc{l — z)] = a{T)n is the density of bosons 
in the condensed (zero-energy) state. Here, a{T) = 
1 - {TlTsEcfl^ where Tbec = 27r[n/C(3/2)]2/3/A/ is 
the critical BEC temperature for non-interacting dilute 
bosons, and C(x) is the Zeta function. In this case, 
n{T) = [Bi/2(z)/A| + Vcnllz]/{n^T), which diverges in 
the thermodynamic limit Il8l | when {N, Vc} oo but 
n — N /Vc is a constant. Thus, we obtain 



T 



EE 



2|pP 
3eo 



Bl/2{Zc) 
B3/2{Zc) 



(n - 



(11) 



which is always smaller than Tbec, and reduces to the 
non-interacting limit of Eq. H1U|) in the absence of BEC. 

The two cases of Bose molecules allows the construc- 
tion of the phase diagrams indicated in Fig.^b) and^^c), 
respectively, where the PBL, FBL, PBEC and FBEC 
are identified depending on the existence of a sponta- 
neous average electric polarization (P) and/or of a BEC 
fraction ((/)). The ferroelectric superfluid phases pro- 
posed here may be experimentally observed with cur- 
rently available cooling techniques only when Tf e is large 
enough. This requirement imposes a condition on the size 
of the electric dipole moments of the molecules, and it is 
discussed next for the Fermi-Fermi Bose molecules. 



Fermi-Fermi Bose molecules: To set the scale, we con- 
sider the specific example of Li-K molecules consisting of 
^Li and "^^K atoms in their ground state, where |p| = 3.6 
Debye |15|. We also choose an equal population mix- 
ture of ^Li and "'"K atoms with parameters N = 10^ 
and V = 10""^ cm^ leading to Tbec ~ 0.099eF and 
Tfe ~ 21Tbec- However, for a molecule with |p| = 1.0 
Debye, Tfe = I-QTbec- Here, e_F = k\'/{2mr) is the 
Fermi energy, where rrir is half of the reduced mass of 
Li and K atoms, and fci? is the Fermi momentum with 
n 
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FIG. 2: Schematic phase diagram of temperature T/Tbec 
versus scattering parameter l/{kpa) for an equal mixture of 
two-species Fermi-Fermi mixtures. 

For the two-species Fermi-Fermi mixtures (such as ^Li 
and '^"K), we propose a schematic phase diagram of tem- 
perature T versus scattering parameter 1/ {kFo) as shown 
m Fig. m In the dia gram, Tq is the pair formation scale 
and Cooper pairs break above this temperature. In addi- 
tion, a two-body bound state does not exist for negative 
a. Thus, a polar bosonic molecule with dipole moment 
IpI ^ may form when both T < Tq and a > criteria 
are satisfied. This also implies that a finite polariza- 
tion Pq may exist below T < Tfe in the phase space 
where T < Tq and a > leading to a possible ferro- 
electric phase. Depending on the relative values of Tfe 
and Tbec, we sketched two cases: (I) Tfe > Tbec and 
(II) Tfe < Tbec which are shown as dotted lines. No- 
tice that as l/{kFa) decreases, the size of the molecules 
increases and p gets smaller, thus producing a reduced 
Tfe- Therefore, an equal mixture of ^Li and ^''K corre- 
spond to case I in Fig. |21 However, other mixtures with 
smaller dipole moments and/or densities may correspond 
to case II. Next, we discuss trapped molecules. 

Heteronuclear molecules in a trap: The trapping po- 
tential of polar molecules is given by ^(r) = — Q;|E(r)p — 
p • E(r), which are due to the coupling of the laser elec- 
tric field E(r) with the induced and permanent dipole 
moments Pind = Q;E(r) and Pper = P, respectively. We 
assume that E(r) = Eox exp[- p'^ / {2wf,) - z^/(2w^)] has 
a gaussian profile controlled by the widths Wp and Wz, 
where = -\- . Thus, the trapping potential can be 
approximated by V{y) = Vq + Amlp^ /2 + Mnlz'^/2, 
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where Vq = —aE^ — |p|£'ocos6', and Qi — {2aE() + 
\p\Eq cos 9)/ (Mwf) are the characteristic frequencies of 
the harmonic trap along i = p,z directions. Here, 9 is the 
angle between p and E(r). When 3> Wp {w^ <C Wp) 
the trap is cigar (disc) shaped. 

In the ground state of the system, the Thomas-Fermi 
approximation can be applied leading to 



Um{v) =p- V{y) + / dr' J(r, r')no(r') 



(12) 



and the electric polarization density is given by P(r) — 
pno(r). The integral equation can be solved analytically 
when (3 — [|pp/(47reo)]/C^ ^ 1, corresponding to a small 
ratio between the characteristic electric dipolar energy 
|pp/(47r£oV'c) and the characteristic contact interaction 
energy U /Vc- To zeroth order in (3, the electric polariza- 
tion is P(r) = PTF(r) = pnTF(r), where 



(13) 



Here, Umax = " Vo)/U, pi = {p ~ VQ)/{MVtl) and 
zl = {p- Vo)/(MVt1). The correction to first order in (3 
for the cigar shaped trap is 

5P(r) = P/3^n^ [ci - C2(pVPc " ^V^c)] , (14) 

where ci and C2 are numerical coefficients, and n = N /Vc 
is the average molecular density. Since E(r) is much 
stronger than the local electric fields produced by elec- 
tric dipolar interactions, the net polarization P(r) = 
PT_F(r) + 5P(r) points along E(r). 



FIG. 3: Direction of the polarization vector during time of 
flight for a ferroelectric state, when the external electric field 
is a) perpendicular, and b) parallel to the easy axis of polar- 
ization. The clouds in each frame are not drawn to scale. 

To distinguish between the fully polarized paraelectric 
superfluid and the true ferroelectric superfluid, we pro- 
pose two experiments, as shown schematically in Fig. |3 
where E(r) is turned off between times I and II. In case 
(a), the external electric field and the easy axis for polar- 
ization are perpendicular. For the ferroelectric state, the 
polarization flips between times I and II (within nanosec- 
onds to microseconds) towards the axial direction which 
corresponds to the easy axis for the electric polarization. 
This flip produces a short pulse of coherent dipolar radia- 
tion proportional to N"^ . However, for the fully polarized 



paraelectric state the polarization relaxes and reduces to 
zero within microseconds, producing a longer pulse of in- 
coherent dipolar radiation proportional to N . A second 
coherent pulse may occur in time of flight for the ferro- 
electric state as the anisotropy inversion in the cloud ex- 
pansion takes place between times HI and IV, and causes 
a second flip of the electric polarization towards the new 
easy axis. In case (b), the external electric field and 
the easy axis for polarization coincide. For the ferroelec- 
tric state, the polarization remains along the easy axis 
between times I and II, and between times HI and IV 
when the anisotropy inversion occurs, the polarization 
flips towards the new easy axis with emission of coherent 
dipolar radiation (proportional to N'^). In this case only 
one pulse (within nanoseconds to microseconds) should 
be observed after a few miliseconds of expansion. For the 
paraelectric state similar results as in case (a) apply. 

In conclusion, we analyzed the possibility of a ferro- 
electric transition in Bose-Bose, Bose-Fermi or Fermi- 
Fermi heteronuclear molecules. This transition is char- 
acterized by the appearance of a spontanous electric po- 
larization below a critical temperature where the dielec- 
tric function diverges. We obtained the order parameter 
equation, evaluated the transition temperature and the 
electric polarization for ultracold heteronuclear (Bose or 
Fermi) molecules. We discussed the existence of a fer- 
roelectric Fermi liquid phase for polar Fermi molecules, 
and the existence of a ferroelectric superfluid phase for 
polar Bose molecules characterized by the coexistence of 
ferroelectric and superfluid orders. We also proposed an 
experiment to detect ferroelectric correlations via the ob- 
servation of coherent dipole radiation pulses during time 
of flight. We thank NSF (DMR-0304380) for support. 
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